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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2938 

ANALYSIS  OF  HEAT  ADDITION  IN  A CONVERGENT -DIVERGENT  NOZZLE 
By  Donald  P.  Hearth  and  Eugene  Perchonok 


SUMMARY 

The  effect  of  heat  addition  to  a divergent  stream  with  initially 
sonic  flow  is  analyzed  and  the  variation  of  exit  Mach  number,  static 
pressure,  and  total  pressure  is  presented.  Application  of  these 
results  to  the  diverging  portion  of  a convergent -divergent  nozzle  indi- 
cated that  nozzle  heat  addition  delays  nozzle  overexpansion  and 
affects  the  jet  thrust  appreciably.  Moreover,  misleading  nozzle  coef- 
ficients will  he  computed  if  heat  addition  in  the  nozzle  is  ignored. 


INTRODUCTION 

It  has  been  shown  (ref.  l)  that  an  increase  in  jet  thrust,  over 
that  achieved  with  a sonic  outlet,  can  be  realized  by  complete  expan- 
sion of  the  exhaust  gases  in  a convergent -divergent  nozzle.  In  the 
calculation  of  this  process  and  in  the  evaluation  of  experimental  noz- 
zle data  the  flow  process  in  the  nozzle  is  usually  assumed  to  be  adia- 
batic. However,  for  nozzles  located  directly  downstream  of  combustion 
chambers  (such  as  ram- jet,  rocket,  or  turbojet  afterburner  nozzles), 
heat  may  be  released  to  the  working  fluid  during  the  expansion  process 
in  the  divergent  portion  of  the  nozzle  as  a result  of  inefficient  or 
delayed  combustion.  The  purpose  of  this  report  is  to  treat  analytically 
heat  addition  to  a divergent  stream  with  initially  sonic  flow  and  to 
employ  the  results  thus  obtained  in  evaluating  the  effect  of  delayed 
combustion  on  convergent -divergent  nozzle  performance.  In  this  analy- 
sis, which  was  conducted  at  the  NACA  Lewis  laboratory,  use  is  made  of 
the  general  equations  developed  in  reference  2 for  heat  addition  to  a 
subsonic  divergent  stream. 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 
A area,  sq  ft 

F jet  thrust,  mV  + A(p  - pQ),  lb 
K constant,  defined  in  equation  (3) 
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M Mach  number 

m mass  flow,  pAV,  slug/sec 

P total  pressure,  lb/sq  ft 

p static  pressure,  lb/sq  ft 

T total  temperature,  °F  abs 

,1* 

V ',^yelocity,  ft/sec 
x constant,,  defined  in  equation  (5) 

X ratio  of  specific  heats 

p static  density,  slug/cu  ft 

1 total-temperature  ratio 

Subscripts : 

c combustion  chamber 

n'  diverging  nozzle 

0 free  stream 

1 intermediate  combustion-chamber  station 

2 combustion>*chamber  outlet 

3 nozzle  exit 

* nozzle  throat 
Superscripts : 

* conditions  for  incremental  heat  addition  prior  to  nozzle 

ANALYSIS 

A schematic  diagram  of  the  case  to  be  considered,  a combustion 
chamber  followed  by  a convergent -divergent  exit  nozzle,  is  shown  in 
figure  1.  Temperature  ratios  of  practical  interest  are  considered  and 
for  analysis  purposes  all  nozzle  heat  addition  is  assumed  to  occur  in 
the  diverging  portion  of  the  nozzle.  That  heat  addition  which  occurs 
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in  the  convergent  portion  of  the  nozzle  is  treated  as  part  of  the 
combustion-chamber  temperature  rise,  thereby  influencing  slightly  the 
total  pressure  available  at  the  nozzle  throat. 


Assumptions 

The  following  assumptions  have  been  made  in  the  analysis : 

1.  The  flow  can  be  considered  as  inviscid  and  one -dimensional. 

2.  The  working  fluid  is  a perfect  gas  with  a ratio  of  specific 
heats  constant  at  a value  of  1.30.  (A  negligible  effect  on  the  com- 
puted results  was  observed  by  using  constant  y values  of  1.32  and 
1.28.  If,  however,  the  value  of  y changes  during  the  process,  the 
results  would  be  somewhat  affected.  For  the  range  of  conditions  con- 
sidered, the  value  of  y increases  during  the  expansion  process.  The 
inaccuracy  introduced  by  the  simplifying  assumption  of  a constant  y 
was  not  considered  of  sufficient  magnitude  to  justify  additional 
refinement  of  the  computations . ) 

3.  The  expansion  process  in  the  convergent  portion  of  the  nozzle 
is  isentropic. 

4.  The  nozzle  throat  is  choked  (M  = 1.0)  at  all  times. 


Heat  Addition  in  the  Divergent  Portion  of  a Choked 
Convergent -Divergent  Nozzle 

An  analysis  of  subsonic  heat  addition  in  a diverging  channel  has 
been  made  (ref.  2),  and,  in  general,  the  equations  developed  therein 
also  apply  to  the  supersonic  case.  In  the  development  of  these  equa- 
tions, the  differential  forms  of  the  equations  for  conservation  of 
momentum  and  energy  are  used  in  conjunction  with  the  equation  of  state 
for  a perfect  gas.  The  following  expressions  relating  the  changes  in 
static  pressure,  area,  total  temperature,  and  Mach  number  can  be 
obtained  (see  ref.  2): 


dp 

P 


yM2 

1-yM2 


^ dM2  \ 


(1) 
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dA  (1  + rM2)dT  (1  - M2)dM2 


(2) 


Simple  integration  of  these  equations  is  possible  if  it  is  assumed 
(as  in  ref.  2)  that 


dA  _ „ dT 
A K T 

The  integrated  form  of  equation  (3)  is: 


(3) 


(3a) 


The  rate  of  heat  addition  may  be  changed  by  altering  the  constant  K, 
the  value  of  -which  is  determined  from  equation  (3a)  by  the  temperature 
rise  considered  and  the  over-all  change  in  flow  area  through  which  the 
temperature  rise  occurs.  Equation  (3)  may  not  necessarily  describe  the 
actual  way  in  which  heat  is  added  but  was  considered  adequate  for 
obtaining  indications  of  trends  in  nozzle  performance. 

When  equations  (l),  ( 2 ),  and  (3)  are  combined  with  the  equation 
for  conservation  of  mass,  the  following  expressions  relating  the  static 
pressure  ratio,  temperature  ratio,  and  Mach  number  change  across  the 
diverging  portion  of  the  nozzle  can  be  derived  (see  ref.  2): 


where 


and 


kl-x 


2K  + rM, 


. 1 - 2K  + fM3' 


r + i + 2K(r-i) 


2 2 1 


(4) 


(5) 


(6) 
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The  ratio  of  total  pressure  at  the  tvo  stations,  * and  3,  is  also 
of  interest  and  is  found  by  expressing  the  total  pressures  in  terms  of 
the  static  pressures  and  Mach  numbers : 


for  the  special  case  "wherein  the  rate  of  area  change  with  heat 
addition  is  such  as  to  maintain  a constant  Mach  number  throughout  the 
heated  region,  the  following  expressions  are  obtained  from  equations  (l) 
and  (2) : 


rM2-! 


J 
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For  this  limiting  case  where  M3  = M^  = 1.0,  the  value  of  K is  1.15 
(f  = 1.30).  For  values  of  K > 1.15,  the  exit  Mach  number,  as  well  as 
the  Mach  number  at  any  point  downstream  of  the  throat,  is  supersonic 
and  greater  than  M*. 


Thrust  Equations 

Jet  thrust  is  defined  as 

F = mV  + A(p  - pQ)  (10) 

If  equation  (10)  is  expressed  in  terms  of  Mach  number  and  total 
pressure,  a ratio  can  be  derived  relating  the  jet  thrust  at  the  outlet 
of  a convergent -divergent  nozzle  to  the  jet  thrust  at  the  nozzle 
throat.  This  latter  value  is  the  thrust  for  a sonic  discharge.  Thus, 


For  M.  = 1.0  and  complete  isentropic  expansion  (P  = PT,  T„  = T~, 
and  P3  = Pq)>  equation  (ll)  reduces  to: 


(12) 


For  this  process,  the  exit  Mach  number  M3  and  the  area  ratio  A-^/a^ 
are  determined  from  the  nozzle  pressure  ratio  P^/pq  b..  making  use  of 
isentropic  flow  relations . 

If,  however,  heat  is  released  to  the  working  fluid  as  it  passes 
through  the  diverging  section  of  the  nozzle,  the  expansion  process  is 
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no  longer  isentropic  and  it  is  necessary  to  solve  equation  (ll)  for 
the  change  in  jet  thrust.  For  given  values  of  x and  A3/A#,  a 
simultaneous  solution  of  equations  (3),  (4),  and  (6)  is  required. 
Appropriate  substitution  in  equations  (7)  and  (ll)  then  yields  the 
desired  thrust  ratio. 

Thus  far,  the  incremental  temperature  rise  has  been  considered  as 
occurring  supersonically  in  a diverging  nozzle.  It  is  also  of  interest 
to  consider  this  incremental  temperature  rise  as  occurring  subsonically 
ahead  of  the  exhaust  nozzle  between  stations  1 and  2.  From  the 
conservation  of  mass,  equation  (13)  can  be  derived  relating  the  con- 
ditions at  the  nozzle  throat  with  and  without  the  incremental 
combustion-chamber  heat  addition  x . 


(13) 


If  equation  (13)  is  used  in  conjunction  with  an  expression  similar 
to  equation  (ll),  a jet  thrust  ratio  for  a sonic  outlet  nozzle  with  and 
without  the  incremental  xc  occurring  prior  to  the  nozzle  can  be 
obtained: 


(i  + r)  - 

* 

(*  • sT 

T 

P,.  / 

t \ r-i 

L(1+r)-^( 

• v) 

(14) 


It  should  be  noted  that  the  absolute  value  of  the  over -all  x at  sta- 
tion 1 before  the  addition  of  the  incremental  ^ c does  not  influence 
the  results  obtained. 


The  solution  of  equation  (14)  for  various  values  of  x0  requires 
thao  the  reduction  in  total  pressure  at  the  nozzle  threat  (p^.1  < P#) 
due  to  the  incremental  subsonic  heat  addition  be  evaluated.  From  the 
conservation  of  mass  and  momentum,  the  following  equations  were  derived 
in  a manner  similar  to  that  presented  in  reference  3: 


Fc 


1 + 


rMl 


4- 


r-i 

2 


T 


(15) 
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M22(l  + g-  Mg2)  Mj2(l  + g:  Ml2) 

(1  + tm22)2  T°  (l  + tm-l2)^ 


(16) 


In  the  examples  computed,  M-^  was  assumed  equal  to  0.5. 

Equation  (14)  indicates  the  gains  in  thrust  at  a sonic  throat  due 
to  the  incremental  heat  addition  made  prior  to  the  nozzle.  If  this 
incremental  temperature  rise  is  accompanied  by  isentropic  expansion  to 
ambient  static  pressure  in  a diverging  nozzle  (eq.  12),  a further 
increase  in  thrust  occurs.  When  this  result  is  compared  with  the 
result  from  equation  (ll)  for  the  same  incremental  t,  the  effect  of 
incomplete  combustion  with  the  resulting  heat  addition  in  the  nozzle 
can  be  determined. 


DISCUSSION  OF  RESULTS 

Heat  Addition  to  a Diverging  Stream  with  Initially 

Sonic  Flow 

Presented  in  figures  2,  3,  and  4 as  a function  of  expansion  ratio 
and  total- temperature  ratio  are  variation  of  exit  Mach  number,  static 
pressure  ratio,  and  total -pressure  ratio  for  heat  addition  to  a diverg- 
ing stream  with  initially  sonic  flow.  The  trends  are  generally  the 
same  as  for  supersonic  heat  addition  in  a constant-area  duct  (ref.  4). 
For  example,  as  shown  in  figure  2,  the  exit  Mach  number  approaches 
unity  as  he'at  is  added  to  the  supersonic  stream.  However,  increasing 
amounts  of  heat  addition  are  required  as  the  expansion  ratio  is 
increased  in  order  to  achieve  an  exit  Mach  number  M3  equal  to  1.0. 

The  effect  of  this  mode  of  heat  addition  on  exit  static  pressure 
is  shown  by  figure  3,  which  indicates  that  the  addition  of  heat  to  a 
diverging  supersonic  stream  causes  an  increase  in  the  exit  static  pres- 
sure. Therefore,  a convergent -divergent  nozzle,  when  operating  at  the 
design  pressure  ratio  for  complete  isentropic  expansion,  becomes  under- 
expanded as  heat  is  released  in  the  divergent  portion  of  the  nozzle. 
Therefore,  to  prevent  overexpansion  when  the  nozzle  is  operating  at 
lower  than  design  pressure  ratios,  heat  may  intentionally  be  added  in 
the  diverging  portion  of  a fixed-geometry  nozzle. 

As  heat  is  added  in  the  manner  being  considered,  a loss  in  total 
pressure  occurs  (fig.  4)  which  is  considerably  greater  than  the  loss 
accompanying  subsonic  heat  addition.  This  loss  in  total  pressure 
increases  as  the  area  ratio  is  increased  because  of  the  succeedingly 
higher  Mach  number  at  which  the  heat  addition  occurs . 
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The  effect  on  jet  thrust  of  expanding  the  flow  from  the  sonic  sta- 
tion and  simultaneously  adding  heat  during  the  expansion  process  is 
shown  in  figure  5.  Presented  (as  determined  by  equation  (ll))  is  the 
ratio  of  the  jet  thrust  at  the  nozzle  exit  to  the  jet  thrust  at  the 
throat  F^/F^,  for  various  temperature  ratios  Tn  and  exit -to -throat 
area  ratios  Aj/A^.  Also  included  on  the  figure  are  lines  indicating 
the  pressure  ratios  P^/p0  required  for  complete  expansion  at  the 
operating  conditions  being  considered. 

The  gain  in  jet  thrust  due  to  expanding  the  flow  isentropically  to 
ambient  static  pressure  as  compared  with  sonic  discharge  is  shown  by 
the  zero  heat  addition  curve,  Tn  = 1.0.  This  is  the  conventional  thrust 
gain  usually  computed  for  convergent -divergent  nozzle  performance  and 
varies  from  approximately  4 percent  at  a pressure  ratio  of  5 to  approxi- 
mately 17  percent  at  a pressure  ratio  of  25.  If,  however,  heat  is 
released  in  the  divergent  portion  of  the  nozzle  (increasing  values  of 
rn),  additional  thrust  gains  result.  At  a of  1.20,  for  example, 

a thrust  gain  from  3 to  4 percent  results  over  the  area  ratio  range 
considered.  It  should  be  emphasized,  however,  that  the  thrust  gains 
indicated  on  this  figure  due  to  increasing  value  of  Tn  occur  only 
because  additional  energy  is  being  added  to  the  working  fluid  and 
that  the  heat  released  prior  to  the  nozzle  remains  unchanged. 

The  following  example  illustrates  the  use  of  figure  5:  A nozzle 

with  a pressure  ratio  P^/pq  = 7*0  and  an  expansion  ratio  A3/a^=  1.70 
(point  A)  is  considered.  At  these  conditions  the  jet  thrust  has  been 
increased  6 percent  over  that  of  a sonic  outlet  by  complete  isentropic 
expansion.  If,  however,  additional  energy  is  added  to  the  working 
fluid  in  the  form  of  heat  released  in  the  diverging  section  of  the 
nozzle  such  that  rn  = 1.20,  the  jet  thrust  can  be  increased  3.5  per- 
cent (point  B)  and  the  nozzle  becomes  underexpanded  (fig.  3).  In  order 
to  achieve  complete  expansion,  the  area  ratio  may  be  increased  to  a 
value  of  2.13  (point  C) . At  this  point  the  jet  thrust  has  been 
increased  13  percent  over  that  of  a sonic  outlet  (a  gain  of  7 percent 
over  the  isentropic  expansion  case)  and  the  nozzle  is  at  its  optimum 
operating  condition  for  the  Tn  and  the  pressure  ratio  considered. 

The  family  of  pressure  ratio  curves  shown  in  figure  5 not  only  indi- 
cates the  area  ratio  variation  required  to  maintain  complete  expansion 
with  nozzle  heat  addition  at  any  given  pressure  ratio  but  also  shows 
the  maximum  thrust  possible  as  compared  with  a sonic  outlet  at  each 
pressure  ratio  if  additional  energy  is  added  to  the  system  by  releasing 
heat  in  the  diverging  portion  of  the  nozzle. 

The  example  shown  in  figure  5 also  illustrates  the  use  of  nozzle 
heat  addition  to  prevent  over expans ion.  A fixed  geometry  nozzle  is 
considered  to  be  designed  for  complete  isentropic  expansion  at  a pres- 
sure ratio  of  approximately  10  and  an  area  ratio  of  2.12  (point  D) . 
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If  the  pressure  ratio  is  reduced  to  7.0  (as  by  a decrease  in  altitude) 
and  since  the  nozzle  area  ratio  is  fixed  at  2.12,  the  nozzle  would  "be 
overexpanded  if  the  flow  were  expanding  isentropically . However,  if 
heat  were  intentionally  added  in  the  nozzle  so  that  “t  = 1.20,  the 
static  pressure  at  the  exit  would  be  increased  to  ambient  static  pres- 
sure so  that  complete  expansion  would  then  exist  (point  C) . 

At  any  given  pressure  ratio,  the  jet  thrust  ratio  increases  essen- 
tially linearly  with  Tn.  Also,  for  a given  value  of  Tn  the  rate  of 
the  increase  in  thrust  ratio  decreases  as  the  nozzle  pressure  ratio  is 
raised.  In  addition,  a larger  area  ratio  change  is  required  to  maintain 
complete  expansion  at  the  higher  pressure  ratios  than  at  the  lower  pres- 
sure ratios . 


Thrust  Comparison  Between  Nozzle  and 
Combustor  Heat  Addition 

A comparison  on  a jet  thrust  basis  between  heat  addition  in  the 
diverging  section  of  a nozzle  and  in  the  combustor  prior  to  the  nozzle 
is  presented  in  figure  6.  The  results  are  presented  in  terms  of  a 
ratio  between  the  jet  thrust  with  incremental  heat  addition  and  the 
jet  thrust  with  a sonic  outlet  and  no  incremental  heat  addition. 

Although  only  results  for  a nozzle  pressure  ratio  Pq/Po  (Pq  = P*  for 
tc  = 1.0)  of  10  are  shown,  similar  trends  were  found  for  all  the  pres- 
sure ratios  considered. 

Case  I in  figure  6 indicates  the  gains  in  thrust  possible  if  the 
additional  energy  were  added  in  the  form  of  heat  released  prior  to  the 
nozzle  and  the  flow  expanded  isentropically  to  ambient  static  pressure 
through  the  nozzle.  If  the  additional  energy  were  added  in  the  diverg- 
ing portion  of  a convergent -divergent  nozzle  designed  for  complete 
isentropic  expansion  (T  = constant),  the  variation  in  thrust  would  be 
as  indicated  b^  case  II.  Since  subsonic  heat  addition  is  more  efficient 
than  supersonic  heat  addition,  the  gains  in  thrust  for  the  same  energy 
input  are  greater  for  case  I than  for  case  II.  However,  this  difference 
can  be  appreciably  reduced  if  the  area  ratios  are  increased  to  maintain 
complete  expansion  as  heat  is  added  in  the  nozzle  (case  III). 

Results  such  as  presented  in  figure  6 have  been  used  to  determine 
the  effect  of  incomplete  combustion  cn  jet  thrust.  In  the  usual  theo- 
retical jet  thrust  calculations  it  is  assumed  that  all  the  energy 
has  been  added  prior  to  the  nozzle  and  that  the  nozzle  expands  the  flow 
isentropically  to  ambient  static  pressure.  Differences  between  measured 
and  theoretical  values  of  thrust  are  then  usually  ascribed  to  inefficien- 
cies in  the  expansion  process  and  are  presented  in  the  form  of  nozzle 
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coefficients.  However,  if  all  of  the  heat  released  by  combustion  is 
not  added  to  the  flow  in  the  combustion  chamber,  but  instead  a portion 
of  it  is  released  in  the  divergent  portion  of  a convergent -diver gent 
exit  nozzle  (as  by  delayed  combustion),  the  usual  nozzle  coefficients 
do  not  adequately  describe  the  efficiency  of  the  nozzle. 

An  example  of  the  difficulty  encountered  in  assigning  coefficients 
to  nozzles  used  with  ram  jets,  afterburners,  and  rockets  is  indicated 
in  figure  7.  The  jet  thrust  was  computed  for  an  ideal  frictionless 
convergent -divergent  exit  nozzle  with  an  area  ratio  corresponding  to 
complete  expansion  with  zero  nozzle  heat  addition  at  each  pressure 
ratio  considered.  It  was  assumed  that  in  one vcase  F^.>  an  incremental 
amount  of  the  total-temperature  rise  across  the  engine,  occurred  in  , 
the  diverging  portion  of  the  nozzle.  In  the  other  case  computed,  F^', 
this  increment  of  heat  addition  was  assumed  to  occur  ahead  of  the 
nozzle  throat.  The  ratio  of  these  two  jet  thrusts  indicates  the 
effect  of  delayed  combustion  and  is  presented  as  a function  of  the 
incremental  total- temperature  ratio  and  nozzle  pressure  ratio  (fig.  7) . 
It  is  apparent  that  if  part  of  the  temperature  rise  occurs  in  the 
nozzle  instead  of  in  the  combustion  chamber,  the  jet  thrust  may  be 
reduced  considerably.  This  thrust  loss  may  be  erroneously  assigned  as 
exit  nozzle  inefficiency  and  the  nozzle  given  misleading  coefficients. 


SUMMARY  OF  RESULTS 

An  analysis  of  heat  addition  to  a diverging  stream  with  initially 
sonic  flow  yields  the  following: 

1.  Heat  addition  in  the  diverging  portion  of  a convergent - 
divergent  nozzle  appreciably  influences  the  jet  thrust  and  must  be 
considered  when  evaluating  nozzle  performance. 

2 . As  heat  is  added  to  a diverging  stream  with  initially  sonic 
fluv,  the  exit  Mach  number  (which  is  supersonic)  approaches  unity,  the 
exit  static  pressure  increases,  the  exit  total  pressure  decreases,  and 
the  total  momentum  at  the  exit  increases. 

3.  Although  heat  addition  in  the  diverging  portion  of  a convergent- 
divergent  nozzle  produces  an  increase  in  jet  thrust,  the  increase  is 
greater  if  the  same  heat  is  added  prior  to  a completely  expanding  nozzle. 

4.  Overexpansion  of  a convergent -divergent  nozzle  may  be  reduced 
by  the  addition  of  heat  in  the  divergent  portion  of  such  a nozzle. 


Lewis  Flight  Propulsion  Laboratory 
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Figure  2.  - Variation  ex:*  Much  n :mt  :r  with  expansion  ratio  for  various  amounts  of  supersonic 
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Figure  6.  - Thrust  comparison  of  several  methods  of  heat  addition 
and  flow  expansion.  P-^/pq  **  10-0. 


Figure  7.  - Effect  of  supersonic  heat  addition  on  Jet  thrust 
of  ideal  nozzle. 
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